Role of Disorder on the Quantum Critical Point of a Model for Heavy Fermions 
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A zero temperature real space renormalization group (RG) approach is used to investigate the role 
of disorder near the quantum critical point (QCP) of a Kondo necklace (XY-KN) model. In the pure 
case this approach yields J c = implying that any coupling J ^ between the local moments and 
the conduction electrons leads to a non-magnetic phase. We also consider an anisotropic version 
of the model (X — KN), for which there is a quantum phase transition at a finite value of the 
ratio between the coupling and the bandwidth, (J/W). Disorder is introduced either in the on- 
site interactions or in the hopping terms. We find that in both cases randomness is irrelevant in 
the X — KN model, i.e., the disorder induced magnetic-non-magnetic quantum phase transition is 
controlled by the same exponents of the pure case. Finally, we show the fixed point distributions 
Pj(J/W) at the atractors of the disordered, non- magnetic phases. 
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Heavy-fermion materials exhibit a variety of interest- 
ing phenomena and different ground states including 
antifcrromagnetic, superconductor, Kondo-insulator and 
metallic []lJ0]- Most of the properties of these systems 
can be attributed to their proximity to a magnetic quan- 
tum critical point (QCP) S. An important and actual 
problem in the study of these materials is to understand 
the role of disorder on the quantum critical point. This is 
also a general problem in statistical mechanics which had 
a lot of progress in the last years mainly due to the study 
of random quantum spin chains 0]. For heavy fermions 
it is particularly relevant to understand the role of dis- 
order since a standard way of bringing these systems to 
the QCP is by doping, which necessarily perturbs the 
system ||. 

The disordered Kondo model || and the Griffiths 
phase approach |j have been proposed to deal with dis- 
order in heavy fermions. In spite of bringing progress 
to this field they are, at this point, mostly phcnomcno- 
logical theories. It is important that results can be ob- 
tained starting from a microscopic Hamiltonian which is 
closest to describe heavy fermions and that can incor- 
porate disorder at a basic level. In this Communication 
we present a non-perturbative real space renormalization 
group approach to the one-dimensional Kondo necklace 
(KN) model at zero temperature, both with and with- 
out disorder. This model was proposed by Doniach [|| to 
study heavy fermions and emphasizes magnetic degrees 
of freedom [P-|l2]| neglecting charge fluctuations. It in- 
corporates the essential physics of these systems, namely 
the competition between Kondo effect and magnetic or- 
dering. The quantum phase transition of this model is 
a magnetic transition which is the actual case in most 
heavy fermions materials. As far as we know the role 
of disorder in the KN model has never been considered 
and the present Communication represents a first step 
towards describing its effects in heavy fermions starting 
at a microscopic level. The KN model is given by the 
Hamiltonian, 
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where and S 11 , fi = x,y,z are spin-1/2 Pauli matrices 
denoting the spin of the conduction electrons and those 
of the local moments, respectively. The sites i and i+1 
are nearest-neighbors on a chain of L sites and Wi is an 
antifcrromagnetic coupling which represents the hopping 
of the conduction electrons between neighboring sites. In 
the present study, either the on-site interactions between 
the conduction electrons and the local moments, Ji > 
or the hopping Wi are uncorrelated quenched random 
variables with initial probability distributions, Pj(Ji/W) 
and Pw(J/Wi) respectively. A closely related model, 
which we also study here is the A-Kondo necklace (X — 
KN) model where the band of conduction electrons is 
represented just by the Ising term, J2i=i Wi(j x ia x i + \. 

In the pure case our zero temperature RG calculations 
yield distinct results for these models. For the X — KN 
we find an unstable fixed point at a finite value of (J/W) 
separating an antiferromagnetic phase from a spin com- 
pensated, Kondo-like phase, which can be fully charac- 
terized. For the XY — KN, any interaction J gives rise 
to a dense Kondo state. The latter result, J c = 0, for the 
XY — KN is in agreement with the most recent density- 
matrix renormalization group p"3| and quantum Monte 
Carlo calculations Q . 

The introduction of disorder in the on-site interactions 
J or in the hopping terms leads us to consider the renor- 
malization of the probability distributions, Pj(J/W) and 
Pw(J/W). For the X — KN model, disorder either in the 
interactions or in the hopping turns out to be an irrele- 
vant perturbation at the pure system fixed point, i.e., the 
zero temperature exponents associated with the transi- 
tion from the antifcrromagnetic (AF) phase to the non- 
magnetic phase in the presence of disorder are those of 
the pure case. This is in strong contrast with the Ising 
model in a transverse field, which has been used as a 
model for heavy fermions (?]], for which disorder is a rel- 
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evant perturbation at d = 1 |L5| and probably also for 
d < 4 Jl6|]. The atractor of the antiferromagnetic phase 
in the presence of disorder is the same of the pure case. 
There is a new atractor associated with the disordered 
non-magnetic phase which is characterized by the fixed 
point probability distributions P}(J/W) and P^(J/W) 
at this atractor in each case. In the XY — KN model 
the flows are always towards non-magnetic phases, which 
are different in the pure and disordered cases. The fixed 
point distributions Pj(J/W) in the non-magnetic phases 
of the X - KN and the XY - KN in the presence of 
disorder are extremely wide. They show huge values of 
(J/W) and negligeable weight for small values of the ra- 
tio (J/W) which are progressively being eliminated. 
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FIG. 1. The cells used in the renormalization group trans- 
formation and which correspond to a length scale factor b — 2. 

In the context of the real-space renormalization group, 
the linear chain may be viewed as b bonds in series; 
b is the scaling factor (6 = 2 in Fig. I). For the 
Hamiltonian given in Eq. 0, the non-commutation as- 
pects are present at the cluster level, in the sense that 
individual spins are not in a definite state. This can 
be dealt with by referring to the density matrix in the 
basis |mim2...m n m' 1 m' 2 ...TO^), where S x \m! i ) = m!/\w!^ 
and a x \m,i) = rrii\mi), and defining the renormalization- 
group (RG) transformation by the mapping of diagonal 
elements only. At zero temperature the density ma- 
trix is essentially the ground state projector. This ap- 
proach has been successfully applied to the transverse 
Ising model |l7) and to quantum spin glasses fls| ]. The 
disorder is incorporated through a statistical renormal- 
ization group (SRG) treatment where one follows the ef- 
fect of a RG transformation on the probability distribu- 
tions of the relevant parameters, instead of forcing them 
to a particular form |jq] . Let us consider the case of 
random interactions J. For a given configuration of on- 
site interactions ({Ji/W}) the RG transformation for the 
d = 1 system (Fig. ) is defined by, 
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where K' = (J/W)' is the renormalized ratio of parame- 
ters in the smallest two-site cell. We have assigned a sin- 
gle index i to each bond in the original cell Ki = Ji/W 
consistent with an uncorrelated distribution for these ra- 
tios. The matrix elements, 

(77117713 1/5' '(K')\mim 3 ) = 



and 
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are obtained by performing the partial trace on the inter- 
nal spins keeping those on the terminal sites (along the 
chain) fixed (see Fig. ). For the renormalized cell, the 
left hand side of Eq (||) provides analytical expressions 
for the primed variables which are given by, 
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For the original, larger cell, the right hand side of 
Eq. (g) is written in a 64 by 64 representation and the cal- 
culations have to be performed numerically. The match- 
ing of these quantities above in the two cells yields in 
principle two RG equations that should allow to treat Jj 
and Wi as independent random variables. It results how- 
ever that these equations reduce to a single one which 
points out that, as concerns the ground state of the sys- 
tem, there is a unique relevant variable in the problem, 
namely, the ratio (J/W). This is clear in the pure sys- 
tem but turns out to be the case also in the presence of 
disorder as is further discussed below 

For convenience let us define, 



Z± = (++ \p\ + +) ± 
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for both b — I and b — 2 cells. We then match the 
quantity Z'_/Z' + = Z-/Z+, which is a scale invariant 
and finally obtain, 



(9) 



which is our basic renormalization group equation. We 
emphasize that although we have considered different 
scale invariants to match in the distinct cells we always 
find the same RG equation or trivial ones . We obtain 
below that the cases of disorder in the interactions or in 
the hopping yield the same phase diagram. These re- 
sults are consistent with the existence of a single random 
variable in the problem, the ratio (J/W), that suffices to 
characterize the ground state phase diagram. Further- 
more, as we take Ki = Ji/W or d = J/Wi as random 
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variables, W or J respectively, are uniform only at the 
level of the cells and not in the chain generated by the 
RG procedure. 

We first consider the pure case without disorder. For 
the X — KN these equations yield a quantum critical 
point at K c = G c = (J/W) c = 0.21 separating an an- 
tifcrromagnctic from a non-magnetic, Kondo singlet or 
spin-liquid phase. The critical exponents associated with 
this transition can be obtained in the usual way. For the 
correlation length, given by £ oc with S = K — K c , 

we obtain the exponent, v = 2.24 ss 9/4. The dynamic 
exponent z is obtained from the finite size scaling re- 
lation for the gaps for excitations in the two cells, i.e., 
(AE/AE') Kc = (L'/L) z = b z , at the quantum criti- 
cal point, where b = 2 is the length scale factor. We 
get z — 1.28. Using the quantum hyperscaling rela- 
tion 2 — a = v(d + z) we can also obtain the exponent 
a = —3.12 which is negative but large in module. The 
non-magnetic phase, for K > K c , is a spin-liquid phase 
characterized by a gap AE = \8\ vz which vanishes at the 
QCP with the gap exponent, s = vz = 2.28. For the 
XY — KN model we find that any finite coupling J gives 
rise to a non-magnetic, Kondo-like phase, i.e., K c = 0. 
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FIG. 2. Schematic phase diagram of the X — KN model 
obtained from our renormalization group method. The ar- 
rows give the flow of the RG equations. The fixed point at 
K c — (J/W) c , a 2 = is that of the pure system. 

The effects of disorder in the on-site interactions in the 
pure Kondo necklace model, are considered taking, for 
convenience, a normalized initial probability distribution, 
P(Ki) = ( J ft'f" 1 //3 a r(a))exp(-/C i //3). This is a single 
peaked distribution with an average value (K) = a/3 and 
a variance a 2 = ((K - (K)) 2 ) = a/3 2 . Notice that K > 0, 
i.e., only positive (antiferromagnetic) values of the on-site 
interactions are considered. We start the iteration of the 
RG equation, choosing three bonds distributed accord- 
ing to P(Ki) above, to feed the recursion relations and 



generate a new value K' . This procedure is repeated up 
to 10 6 times and yields a new renormalized distribution 
P'(K') which will be used to feedback the recursion rela- 
tion in the next step of the renormalization process. The 
evolution of the distribution is then followed as the renor- 
malization process goes on. We emphasize that the exact 
form of the initial probability distribution is not relevant 
and the different starting distributions we used always 
evolve to the same fixed point form associated with the 
atractor of the corresponding phase. 

We next consider adding randomness to the on-site 
interactions of the pure X — KN model. The zero- 
temperature RG flow diagram as a function of (K) = K 
and a 2 is shown in Fig. 2. In the presence of disorder the 
system still exhibits two phases: the Ising antiferromag- 
netic phase (X — AF) and a disordered, non-magnetic 
phase. The atractor of the antiferromagnetic, X — AF, 
phase in the presence of disorder is the same of the pure 
AF phase. The probability distribution in this ordered 
phase is single peaked and becomes sharper, i.e., a 2 — > 
as the mean value approaches the value (K) = un- 
der successive renormalizations (see Fig. 3). On the other 
hand, it is clear from the phase diagram in Fig. 2 that the 
atractor of the non-magnetic phase in the presence of dis- 
order is different from the atractor of the non-magnetic, 
spin-liquid phase of the pure case. 
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FIG. 3. The probability distributions Pj(J/W) at two su- 
cessive steps of the statistical renormalization procedure (/ 
and II respectively) as the system iterates to the atractor of 
the disordered, non-magnetic phase (a) and to the atractor of 
the antiferromagnetic phase (6) of the X — KN model. The 
insets show the initial points of the RG processes. The critical 
boundary in Fig. separates the two types of flow shown here. 

The phase transition to the ordered magnetic phase in 
the presence of disorder is governed by the same critical 
exponents of the pure case that we have calculated before. 
This is a consequence that the flow along the critical line 
is towards the semi-stable fixed point of the pure system 
at K c — 0.21, a = which therefore determines the uni- 
versality class of the transition, even in the presence of 
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disorder. The irrelevance of disorder in the renormaliza- 
tion group sense is consistent, according to Harris crite- 
rion, with the negative value of the exponent a = —3.12 
we have obtained previously. The critical line close to 
the QCP vanishes as, a 2 oc \5\ = \K — K c \. The ex- 
act form of the critical frontier in the region represented 
by the dashed line is not numerically accessible. Most 
probably it approaches asymptotically the axis K = 
as there should be no critical value for this type of dis- 
order to destroy the X — AF phase at K = 0. Notice 
that disorder is detrimental to antifcrromagnetism, re- 
ducing the region of the phase diagram where this can 
be found. We have also carried out the statistical RG 
procedure for the case of disorder in the hopping terms 
using the same initial distributions, but now for the ran- 
dom variable Gi — J/Wi. It turns out that the same 
phase diagram of Fig. 2 is obtained. 

The effects of disorder in the on-site interactions on 
the XY — KN, which in the pure case has only a non- 
magnetic Kondo phase, are quite distinct, from the RG 
point of view, compared with those in the anisotropic X— 
KN model. Firstly, disorder is a relevant perturbation 
at the K c = fixed point since for any initial values 
of (K) ^ and a 2 , the renormalized variance always 
increases in the renormalization process. Next, the flow 
of the distribution P(J/W) towards the atractor of the 
random Kondo phase is much faster than for the X — 
KN. The phase diagram of the XY — KN has only 
two atractors, that of the pure Kondo phase at a = 0, 
J/W = co and that of the random Kondo phase at a — 
co, (K) — co. This is the first time this result has been 
found using a real space RG. The same considerations 
apply for the case of off-diagonal disorder, i.e., when the 
hoppings are random, with the replacement (K) — > (G). 

In the non-magnetic phases, of both the disordered 
X — KN and XY — KN models, for the case of dis- 
order in the interactions Ji, the distributions Pj(K) be- 
come extremely wide under iteration as shown in Fig. 3. 
Small values of K are being eliminated and huge values 
for this ratio are being generated. Note that this behav- 
ior of Pj(K) does not imply the absence of low energy 
excitations since large values of K can occur for J finite 
and W — > jll). In fact the only meaningful variable 
in the present approach is the ratio of these quantities. 
These disordered phases are better characterized by the 
distribution of gaps for excitations which however can 
not be obtained in the present approach since it does not 
take into account excited states. 

In summary, we have studied the effect of disorder 
in the KN model for heavy fermion systems. The 
anisotropic version of the model, the X — KN, is appro- 
priate to describe systems where the ordered magnetic 
phase has a strong Ising character. In this case for pure 
systems we have found a quantum critical point, at a 
finite value of (J/W), separating the ordered magnetic 
phase from the spin-liquid. Remarkably, we have ob- 
tained that at d — 1, where the effects of disorder are 
stronger, this turns out to be an irrelevant perturbation 



at the pure system quantum critical point. This holds 
for both cases of initial disorder in the interactions Ji 
or in the hoppings Wj. The method used here has the 
great advantage of being non-perturbative which allows 
to the treat the effects of weak disorder in the neigh- 
borhood of the quantum critical point of the pure sys- 
tem. Our results are in agreement with the robustness of 
the quantum critical point approach description of non- 
Fermi liquid behavior in systems where disorder is clearly 
present. The present RG method characterizes the dif- 
ferent phases in terms of their atractors. In particular, in 
the disordered phases the same atractor is reached inde- 
pendently of the point of the phase diagram we start the 
statistical iteration procedure. Consequently the present 
approach does not yield information on the existence of 
rare clusters of the incipient magnetically ordered phase 
which characterizes the so-called Griffiths phases. This is 
a feature of the method which considers only the ground 
state of the cells and no excited states are involved in 
the calculations. Another consequence of this feature is 
the existence of a single variable, the ratio (J/W) which 
suffices to characterize the phase diagram of the system 
even in the presence of disorder. This approach however 
can be extended to include excited states and to calcu- 
late the distribution of gaps in which case we expect to 
find evidence of Griffiths phases M] . 
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